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Abstract 

In this paper, we provide a new and effective approach for studying super-exponential 
solution of a retrial supermarket model with Poisson arrivals, exponential service times 
and exponential retrial times and with two different probing-server numbers. We de- 
scribe the retrial supermarket model as a system of differential equations by means 
of density-dependent jump Markov processes, and obtain an iterative algorithm for 
computing the fixed point of the system of differential equations. Based on the fixed 
point, we analyze the expected sojourn time that a tagged arriving customer spends 
in this system, and use numerical examples to indicate different influence of the two 
probing-server numbers on system performance including the fixed point and the ex- 
pected sojourn time. Furthermore, we analyze exponential convergence of the current 
location of the retrial supermarket model to the fixed point, and apply the Kurtz 
Theorem to study density-dependent jump Markov process given in the retrial super- 
market model, which leads to a Lipschitz condition under which the fraction measure 
of the retrial supermarket model weakly converges to the system of differential equa- 
tions. This paper arrives at a new understanding of how the workload probing can 
help in load balancing jobs in retrial supermarket models. 

Keywords: Randomized load balancing, retrial queue; supermarket model, super- 
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exponential solution, density-dependent jump Markov process, fixed point, sojourn 
time, exponential convergence, Lipschitz condition. 

1 Introduction 

This paper analyzes a retrial supermarket model with Poisson arrivals, exponential service 
times and exponential retrial times and with two different probing-server numbers. It 
proposes a new analytic approach for discussing retrial queueing networks. This new 
approach is based on both differential equation approximation for Markov chains and 
limiting structure for population processes, e.g., see Kurtz |16j . Wormald [43J and Darling 
and Norris [6]. A key result of this paper is an iterative algorithm for computing the 
fixed point, which is always important in the study of supermarket models. Based on the 
fixed point, this paper uses numerical examples to indicate different influence of the two 
probing-server numbers on system performance including the fixed point and the expected 
sojourn time. 

Recently, a number of companies (e.g., Amazon, Google etc.) are offering cloud com- 
puting service and cloud manufacturing technology. This motivates us to study efficient 
load balancing in large-scale networks with many computational and manufacturing re- 
sources. It is well-known that optimal task assignment is a key design issue when one can 
use different resources of large-scale networks. For such task assignment, there are two 
main approaches: Push strategies and pull strategies. Push strategies are mainly used in 
centralized systems, such as data centers allocate jobs to different processors; while pull 
strategies are applied to send work from an overloaded processor to underloaded proces- 
sors, or an underloaded resource will ask for work from overloaded resources. In the study 
of push strategies, randomized load balancing is an important method, and specifically, 
retrial supermarket models provide a crucial task assignment mechanism in randomized 
load balancing schemes. 

Randomized load balancing, where a job is assigned to a server from a small sub- 
set of randomly chosen servers, is very simple to implement. It has surprisingly good 
performance (for example, reducing collisions, waiting times, backlogs) in a number of 
applications, such as, data centers, hash tables, distributed memory machines, path selec- 
tion in networks, and task assignment at web servers. Supermarket models are extensively 
used to study randomized load balancing schemes. In the study of supermarket models, 
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a key result by Vvedenskaya, Dobrushin and Karpelevich [39] indicated that when each 
Poisson arriving job is assigned to the shortest one of d > 2 randomly chosen queues with 
exponential service times, the equilibrium queue length can decay doubly exponentially 

in the limit as the population size n — > oo. Also, the stationary fraction of queues with 

d k -i 

at least k customers is p d - 1 . This is a substantially performance improvement over the 
case for d = 1, where the tail of stationary queue length distribution in the corresponding 
M/M/l queue is p k . Further, such an exponential improvement is illustrated by Li and 
Lui [19] for studying doubly exponential solution of the supermarket model with Marko- 
vian arrival processes and phase-type service times, and by Luczak and McDiarmid [22] 
for analyzing the maximum queue length in the supermarket model with Poisson arrivals 
and exponential service times. 

In the past ten years, supermarket models have been formulated by queueing theory 
as well as Markov processes. Much recurrent research dealt with the supermarket model 
with Poisson arrival processes and exponential service times by means of density-dependent 
jump Markov processes, which are applied to describe limiting behavior of the supermar- 
ket model as an infinite-size system of differential equations whose fixed point is doubly 
exponential. Readers may refer to Vvedenskaya, Dobrushin and Karpelevich [39J and 
Mitzenmacher [28]. Certain generalization of supermarket models has been explored in, 
for example, studying simple variations by Vvedenskaya and Suhov |40j . Mitzenmacher and 
Vocking [35], Mitzenmacher [29j [30], [33] , Azar, Broder, Karlin and Upfal [3], Vocking [38] 
and Mitzenmacher, Richa, and Sitaraman |34j : considering non-Poisson arrivals or/and 
non-exponential service times by Li, Lui and Wang |20|. 121]. Li and Lui [19] . Bramson, Lu 
and Prabhakar [3] and Li [18] ; discussing load information by Mirchandaney, Towsley, and 
Stankovic [36], Dahlin [5] and Mitzenmacher |32} 134] : mathematical analysis by Graham 
[11] fl~2] [T3] , Luczak and Norris |24] and Luczak and McDiarmid [22 \ [23] ; using fast Jackson 
networks by Martin and Suhov [26], Martin [25] and Suhov and Vvedenskaya [37] . 

Retrial queues are an important mathematical model for studying telephone switch 
systems, digital cellular mobile networks, computer networks and so on. During the last 
two decades, considerable attention has been paid to the study of retrial queues, for 
example, by survey papers of Yang and Templeton [33], Falin [8], Kulkarni and Liang 
[15] . Artalejo [1] and Gomez-Corral |10j . and by three books of Falin and Templeton [9], 
Artalejo and Gomez-Corral [2] and Li |17j . Note that the study of retrial queueing networks 
is difficult and challenging with few available results up to now, this paper proposes a novel 
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approach for studying a class of important retrial queueing networks. 

The main contributions of the paper are threefold. The first one is to provide a novel 
approach to study a retrial supermarket model with Poisson arrivals, exponential service 
times and exponential retrial times and with two different probing-server numbers. We 
describe the retrial supermarket model as a system of differential equations by means of 
density-dependent jump Markov processes, and give an iterative algorithm for computing 
the fixed point which is shown to be super-exponential. The second one is to give an 
iterative algorithm for computing the expected sojourn time by means of the fixed point, 
and use numerical examples to illustrate different influence of the two probing-server 
numbers on system performance including the expected sojourn time. The third one is in 
analyzing exponential convergence of the current location of the retrial supermarket model 
to the fixed point. Not only does the exponential convergence indicate existence of the 
fixed point, but it also explains that such a convergent process is very fast. Furthermore, 
to study limiting behavior of the retrial supermarket model as its population size goes to 
infinity, we apply the Kurtz Theorem to study density-dependent jump Markov process 
corresponding to the retrial supermarket model. This leads to a Lipschitz condition under 
which the fraction measure of the retrial supermarket model weakly converges to the 
system of differential equations. Different from previous works in the literature of retrial 
queues, this paper gives a more general analytical approach in the study of retrial queueing 
networks. We believe our approach can open a new door to performance evaluation of 
more complicated retrial queueing networks. 

The remainder of this paper is organized as follows. In Section 2, we describe a re- 
trial supermarket model with Poisson arrivals, exponential service times and exponential 
retrial times and with two different probing-server numbers, and derive a system of dif- 
ferential equations based on density-dependent jump Markov processes. In Section 3, we 
provide an iterative algorithm for computing the fixed point of the system of differential 
equations, and show that the fixed point is super-exponential. In Section 4, we provide an 
effective algorithm for computing the expected sojourn time, and use numerical examples 
to indicate different influence of the two probing-server numbers on system performance 
including the expected sojourn time. In Section 5, we study exponential convergence of 
the current location of the retrial supermarket model to the fixed point. In Section 6, we 
apply the Kurtz Theorem to study density-dependent jump Markov process correspond- 
ing to the retrial supermarket model, which leads to a Lipschitz condition under which 
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the fraction measure of the retrial supermarket model weakly converges to the system of 
differential equations. Some concluding remarks are given in Section 7. 

2 Retrial Supermarket Model Description 

In this section, we describe a retrial supermarket model with Poisson arrivals, exponential 
service times and exponential retrial times and with two different probing-server numbers, 
and derive a system of differential equations based on density-dependent jump Markov 
processes. 

Let us describe the retrial supermarket model. Customers arrive at a queueing system 
of n > 1 servers as a Poisson process with arrival rate nX for A > 0, and the service 
times at each server are exponential with service rate [i > 0. There is no waiting space 
before each server, but there is a total orbit of infinite size for all waiting customers. 
Each primary arriving customer chooses d\ > 1 servers independently and uniformly at 
random from the n servers. If there is no customer in one of the d\ chosen servers, then it 
enters the idle server immediately and receives service; otherwise it enters the orbit and 
makes a retrial at a later time. Each retrial arriving customer chooses d,2 > 1 servers 
independently and uniformly at random from the n servers, if there is no customer in one 
of the c?2 chosen servers, then it enters the idle server immediately and receives service; 
otherwise it enters the orbit again and makes a retrial at a later time. Returning customers 
behave independently of each other and are persistent in the sense that they keep making 
retrials until they receive their requested service. Successive inter-retrial times of each 
customer in the orbit are exponential with retrial rate nd for 9 > 0. We assume that all 
the random variables defined above are independent of each other, and that this queueing 
system is operating in the region p = X/fj, < 1. 

In the retrial supermarket model, Figure 1 provides physical structure of the retrial 
supermarket model, and Figure 2 shows that each primary arriving customer and each 
retrial arriving customer can choose d\ and di servers independently and uniformly at 
random from the n servers, respectively. We call d\ and c?2 the two probing-server numbers. 
Note that the two probing-server numbers have different influence on system performance 
including the fixed point and the expected sojourn time. 

The following lemma, which is stated without proof, intuitively provides a sufficient 
condition under which the retrial supermarket model is stable. Note that this proof can be 
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No waiting room before each server 
Figure 1: Physical structure of the retrial supermarket model 




Figure 2: The primary and retrial customers probe the loading of di,d,2 > 1 servers, 
respectively 
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given by a simple comparison argument with the queueing system in which each customer 
queues at a random server (i.e., where d = 1). When d = 1, each server acts like an M/M/l 
retrial queue which is stable if p = A/ p < 1, see Falin and Templeton [9] or Artalejo and 
Gomez-Corral [2]. The comparison argument is similar to those in Winston [42J and 
Weber [H] , thus we obtain two useful results: (1) the shortest queue is optimal due to the 
assumptions on Poisson arrivals, exponential service times and exponential retrial times; 
and (2) the size of the longest queue in the retrial supermarket model is stochastically 
dominated by the size of the longest queue in a set of n independent M/M/l retrial queues. 

Lemma 1 For the retrial supermarket model with Poisson arrivals, exponential service 
times and exponential retrial times and with two different probing- server numbers, it is 
stable if p = X/fi < 1. 

For k > 0, we denote by it) the numbers of busy servers with at least k customers 
in the orbit at time t > 0, and ni (t) the numbers of idle servers with at least k customers 
in the orbit at time t > 0. Clearly, (t) + Uq (t) = n and < (i) , n£ (t) < n 
for k > 0. 



For k > 0, we write 



(W) ( U +\ - n k (*) 



n 



and 



n 

which are the fractions of busy servers and of idle servers with at least k customers in the 
orbit at time t > 0. Let 

X n (Q,t)=xW (0,i) 

and for A; > 1 

In(M) = (xr' (A;- 1,0,4° 
X„(i) = (X n (0,t),X n (l,t),X n (2,t),...). 

The state of the retrial supermarket model can be described as the vector X n (t) for t > 0. 
Since the arrival process to this queueing system is Poisson and the service or /and retrial 
times of each customer are all exponential, {X n (t) ,t > 0} is a Markov process whose state 
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space is given by 

={ (<?i 0) , gP , 9^ ■■): g^ e (0, 1) , > g^ > 0, and ng^ is 
a two-dimensional row vector of nonnegative integers for k > 1}. 



For k > 0, we write 



and 



Let 



4 W) (n,t) = E[xW (k,t) 
5 (n,i) = 4 (M) 



and for k > 1 



Sk (n, t) = (s^l (n, t) , 4 7) (n, t) 

S (re, t) = (S (n, t) , Si (n, t) , S 2 (n, t) , . . .) . 

As shown in Martin and Suhov [26] and Luczak and McDiarmid [22J, the Markov process 
{X n (t),t>0} is asymptotically deterministic as the population size n — > oo. Thus, 



lim„ 



E 



Xn ^ (k, i) 



and limn^ooE 1 



>(k,t) 



always exist by means of the law of 



large numbers for k > 0. Based on this, for k > we write 



and 



S ( k I) (t)= lim ^(n.t), 
n— s>oo 

S k (t) = lim S'jfc (re, t) 



S(t) = (S (t),S 1 (t),S 2 (t). 



Let X (t) = lim 

n— s-oo X n (t). Then it is easy to see from the Poisson arrivals, the exponential 
service times and the exponential retrial times that {X (t) ,t > 0} is also a Markov process 
whose state space is given by 



n 



{ (glWKgW, ...): € (0, l),gW > g^ > for k > l} . 



If the initial distribution of the Markov process {X n (t) ,t > 0} approaches the Dirac delta- 
measure concentrated at a point cj £ then X (t) = lim^—^Q^ X n (t) is concentrated on the 
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trajectory S g = {S (t) : t > 0}. This indicates a law of large numbers for the time evolution 
of the fractions of busy servers and of idle servers. Furthermore, the Markov process 
{X n (t) ,t> 0} converges weakly to the fraction vector S (t) = (So (t) , S\ (t) , S2 (t) , . . .) 
as the population size n — > 00, or for a sufficiently small e > 0, 

lim P{\\X n (t)-S(t)\\>e} = 0, 

n— >oo 

where ||a|| is the Loo-norm of vector a. 

The following proposition shows that the three sequences js'^'' (t) + S^ , {s^"* 
and I S^ (£) j are monotone decreasing, while its proof is easy by means of the definition 
of the fraction measure S (t). Note that the crucial number Kd lt d 2 wu l be indicated in 
Tables 1 to 3 later. 

Proposition 1 For k > 1 

sS3 (*) + 4% a) < s r ^ + ( t) < s r ^ + ( t) = 1. 

At the same time, there exists an integer Kd 1; d 2 > such that for k > Ka lt d 2 
and 



s { k I l 1 ( t )< s k'(t)- 

Now, we set up a system of differential equations in order to determine the fraction 
measure S (t) by means of density-dependent jump Markov processes. To that end, we 
use Figure 3 to provide a concrete example for k > 1 to indicate how to derive the system 
of differential equations. 
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Figure 3: State transitions of a retrial queue in the retrial supermarket model 
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Consider the retrial supermarket model with n servers, and determine the expected 
change in the number of busy servers over a small time period of length dt. Dur- 
ing the time interval [0,dt), the probability that any primary arriving customer joins 

r (w) - • i 

the orbit with k — 1 customers is given by nX ■ 



CM) 



dt, the probability that 



any primary arriving customer enters an idle server with k customers in the orbit is 



given by nX 



AW) 



n,t) 



4 (n,t) /■dt, the probability that any retrial ar- 



riving customer enters an idle server with k + 1 customers in the orbit is given by 



(k + 1) n6 4+i( n >^) dt, and the probability that any customer finishes its required 



d 2 



service and leaves this system with k customers in the orbit is given by n 
Therefore, we obtain 



( w ) I +\ 
fis K k (n,t) 



dt. 



dE 



n 



(W) 



(t) 



dt 



which leads to 



-nX 



■23 <m) 



+ nX 



Si) 



(n,t) 



nX 



+ n 



4 W) (n,t) ' + (k + l)nd 



-fis^ (n, t) 



d 2 



d*]T(M) 
dt 



•22 (m> 



di 



4 W) (^*)J " + + 1)' 
fisf } (n, t) . 



<h 



+ X 

SW) 



Si) 



(n,t) 



'fc+i 



(n,t) 



<i 2 



Noting that linin^oo si'"' ; (n, t) exists for k > and taking n 



(1) 

co in both sides of 



Equation ([I]), we have 

MP (t) 
dt 



s£5(*) 



+ A 



sf*>(t)] l + (A; + l)0[s«(i) 



d 2 



(2) 



Using a similar analysis to Equation ([I]) , we obtain the system of differential equations 
satisfied by the fraction vector S (n, t) = (So (n, t) , S\ (n, t) , S2 (n, t) , . . .) as follows: 



d4 J) (n,t) 



d4 W) (M) 



dt 

for > 1 
ds 



SI) 



dt 



(n,t) 



-X 



Si) 



%'(n,t) '+^s W) (n,t), 



di 



SW) 



(W) 



n,t) 



dt 



SW) 



+ A 



>(n,t) 

4 7) (M) 



di 



(n,t) 



4 W) (M)J ' + (k + l) 



+ 



4 7) ( n > *) - ^4 W) ( n > *) 



<2 2 



(3) 
(4) 



di 



Si) 
'fc+l 



(n,t) 



d 2 



(5) 
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and 



d4 7) (n,t) 
dt 

Si) r„ + \ , .(wi 



-A 



4 7) M) 



£;0 



4 7) (n,t) ■ + A xa^(n,t). 



d 2 



(WO 



(6) 



Note that s K / (n,t) + ' (n,t) = 1. 

Since lim n _ i . 00 (n, t) exists for k > 0, taking n — )• oo in both sides of the system of 
differential vector equations ([3j) to ([6]), we obtain a system of differential vector equations 
for the fraction vector S (t) = (So (t) , Si (t) , S2 (t) , ■ ■ ■) as follows: 



~dJ 



for k > 1 



~dt~ 



dSl w >(t) 



asP (t) 

dt 



423 (*) 



-A 



s^ (t)] 1 + (t) , 



sP (t) 



sf ] (t) - + (k + i)e sl<> (t) 



+ 



di 



da 



<i2 



+ A 



WW -vsrv 



(7) 
(8) 



0) 



and 



dt 



A;0 



Note that sJ J) (t) + (*) = 1- 



(10) 



3 Super-Exponential Solution 

In this section, we provide an iterative algorithm for computing the fixed point of the sys- 
tem of differential equations ([7]) to (|1U|) . and show that the fixed point is super-exponential. 

Let 7To = 7Tq^ and ir^ = (n^l i 71 "^) ^ or k > 1. A row vector tt = (tto,tti,tt2, . . .) is 
called a fixed point of the fraction vector S (t) if lim^ +00 S (t) = tt. It is obvious that if 
7T is a fixed point of the fraction vector S (t), then 



lim 

t— s-+oo 



0. 



Taking t — > +00 in both sides of the system of differential equations (JT]) to (|1Q|) . 
together with S^ (t) + S Q W) (t) = 1, we obtain a system of nonlinear equations as follows: 



vr (/) + tt {W) 

TTq -f- 7T 



I- 



(11) 



11 



7T, 



(J) 



, (WO n 



7T, 



CO 



di 



7T, 



(WO 



+ 



7T 



(0 



d 2 



^6 



0. 



for k > 1 



7T 



(WO 
fc-i 



7T 



(WO 



rfl 



+ (& + !) 



7T 



(0 ' 
fc+1 



d-2 



+ A 



7T 



(/) 



di 



^k 



(W) 



and 



7T 



(/) 



fc0 



7T 



U) 



0. 



(12) 
(13) 

(14) 
(15) 



The following proposition shows that the three sequences j ttJ^ + vr^ | , 1 7r£ V ^' ) | and 
i 71 "!^} are monotone decreasing, while the proof is easy in terms of Proposition [TJ Note 
that the crucial number K^d^ win be indicated in Tables 1 to 3 later. 



Proposition 2 For k > 1 



.(wo , m ^ (wo . m ^ rwo . rr> 



+ vr fc+1 < vr fc + vr fe < vr + tt, 







vlt i/ie same time, i/iere exists an integer K^d^ — swc/i i/iai /or A; > 

(W) JW) 

n k+l 71 k 



and 



7T (/) <7T (/) 



By using the system of nonlinear equations (jlip to (|15p . we begin to derive the closed- 
form solution for the fixed point tt = (ttq, ttx,tt2, •••)■ To that end, we denote by 77* = 1 — 77^ 
a positive solution to the nonlinear equation 



Xx dl + j2X - jj, = 



(16) 



or 



sc = 1 — 1 . 



Proposition 3 There exists the unique positive solution 77* to the nonlinear equation 
Xx dl + fix — [1 = for x £ (0, +00), and 77* € (0, 1). 
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Proof: If rf is an any positive solution to the equation x = 1 — px dl for x G (0, +00), 
then it is clear that 77* G (0, 1) due to the fact that x = 1 — px dl < 1. Let G (x) = px dl + x 
and H (x) = 1. Since G (0) =0 and G(l) = 1 + p > 1, there the two continuous 
functions: y = G (x) and y = 1, must have an intersection point rf G (0, 1) by means of 
the intermediate value theorem, that is, the equation x = 1 — px dl must exist a positive 
solution in the interval (0, 1). Note that ^G(x) = pdix^ 1 + 1 > for x G (0, 1), thus 
is strictly monotone increasing for x G (0, 1). This shows that the positive solution 
rf is unique for the nonlinear equation G (x) = H (x), or x = 1 — px dl . This completes 
the proof. ■ 

Remark 1 The numerical implementations are used to illustrate that the unique positive 
solution rf to the nonlinear equation x = 1 — px dl for x G (0, 1) has two monotone 
properties: rj* is a monotone increasing function of d\ > 1 for a given p G (0, 1), while it 
is a monotone decreasing function of p G (0, 1) for a given d\>\. 

When d\ = 3, the positive solution is 

p= 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

r,* = 0.92169 0.86883 0.82905 0.79728 0.77091 0.74844 0.72890 0.71165 0.69624 

When d\ = 5, the positive solution is 

p= 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

T]* = 0.93031 0.88896 0.85937 0.83633 0.81747 0.80151 0.78770 0.77554 0.76468 

When d\ = 10, the positive solution is 

p= 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

n* = 0.94387 0.91643 0.89785 0.88375 0.87236 0.86281 0.85457 0.84734 0.84089 

When d\ = 20, the positive solution is 

p= 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

r]* = 0.95779 0.94088 0.92989 0.92169 0.91513 0.90967 0.90497 0.90086 0.89719 

When d\ =50, the positive solution is 

p= 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

7]* = 0.97366 0.96547 0.96033 0.95656 0.95357 0.95109 0.94897 0.94711 0.94547 

Let 

<5 = p{l-r] dl ) dl 
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and for k > 1 



The following theorem provides an iterative algorithm for computing the fixed point 
which is shown to be super-exponential. 

Theorem 1 The fixed point it = (7ro,7ri,7T2, ...) is given by 

7r o /) 



7T, 



(W) 



<5o, 



/or > 1 



and 



TT 



A rdi ^ J 



kO 



'jfe-1 



TT 



Proof: By induction, we can derive the above desired result. 
It follows from Equations ([Til . (fT2j) and ([16]) that 



A 



7I"i 



di 



+ /i7r^ / ' 1 — = 0. 



This gives 



and 



(I) * 
TT y ' =i] =1- Vdl 



It follows from (|13p that 

A(l- % ) dl -A^(i-%) d? + 
which follows 



7T 



(J) 



c/2 



fj,p{l-r] dl ) 1 = 



7T 



Together with (|15p for A; = 1, this gives 



TT 



(W) 



pp^ (1 - Vdl f + p 



-5 dl 



■'2 



(17) 



ft. 



14 



Now, we first assume that Equations (|Mj) and ([15]) are correct for the case with I = k, 
that is, 



7T 



CO 



A sdi ^ " 



k6 



fc-1 



and 



Then we need to check that Equations (|14p and (|15p are also correct for the case with 
I = k + 1. To that end, it follows from 03J) that for Z = & + 1 



Ati"A^ + (Hl) 



7T 



(o 

fc+1 



d 2 



A 



d 2 



+ A ( ^C-i ) - ^ = o- 



Noting that 



we obtain 



hence we have 



Mfc = A^i + A ( 



7T 



(0 
fc+1 



A 



(A; + 1) 6 



which, together with (I15p . follows that 



7T 



fc+1 



7T 



(0 
fc+1 



+ P 



7T 



fc+1 



di 



(fc+l)0 



(A; + 1)0 fc 



Mfc 1 + P 
^fc+i- 



A 



(A; + 1) 9 



This completes the proof. 

The following corollary sets up a useful relation between 7r^, and 7rj^ for k > 1. 



Corollary 2 For > 1 



7T 



(0 



d 2 



7T 



(WO 
fc-i 



and there exists an integer K > 1 suc/i i/iai /or k > K 



7T 



(/) 



(fa 



< 



7T 



fc-1 



15 



Furthermore, we have 



lim 

k— >oo 



TT 



(J) 



TT 



(W) 
k-1 



-\di 



0. 



Proof: It follows from (fT7|) that for k > 1 



7T 



<fe _ A 



7T 



(WO 
fc-i 



which leads to 



and 



k6 



TT 



(I) 



TT 



(WO 
fc-i 



7T 



A 



7T 



(W) 



fc-1 



*1 



Hence we have 



lim 

k— ¥00 



TT 



(I) 



d 2 



TT 



(W) 
k-1 



lim 



A 



dl fc-s>oo k8 



TT 



(I) 



d 2 



0. 



< 



AW) 



'k-i 



for k > K due to 



Taking an integer K > X/9, we obtain that 
X/k9 < 1. This completes the proof. ■ 
Using Theorem (H we now consider a special case with d\ = di = d > 1. The following 
corollary provides an explicit expression for the fixed point with d\ = d<i = d. 



Corollary 3 If di = d% = d > 1, then the fixed point tt = (ttq, 7Ti, tt2, •••) is given by 

CO 1 



=P{ 1 -Vd) ■ 



for k > 1 

(W) _ A + k6 



TT, 



k0 



A + (fc-1). 
(fc- 1)0 



and 



7T 



(J) _ / a y A + (fc-i). 



fc^y (A; - 1) i 



A + (fc - 2) 
{k — 2)9 



\ + (k-2)6 
(k -2)6 



d 2 



A + 



p d "l (1 - 7/ d ) 



(18) 



P^il-mf . (19) 
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Proof: By induction, we can easily derive the above desired result. 
If d\ = d.2 = d > 1 , then 



6 = p(l- n d ) c 



and for k > 1 



-p- 



x + ke 
i c 

kO 
X + k9 
k6 



k-l 



A + (fc-l) 
(k-l) 6 



°k-2 



x + ke 
kl 



x + (k-i)e 

(k - 1)6 



A+(fc-2)i 
(fc - 2) 6> 



id 2 



A + . 



This gives 



7T 



A + fc0 



A + (fc- 1). 
(fc- 1)0 



A + (fc - 2) 
(k-2)6 



d 2 



and 



fc^) " 

A\ h + (fc-i) 
fc^y (fc - 1)0 

This completes the proof. 



A + (fc-2)fl 
(k-2)9 



x + e 



p d - 1 {i-m) 



A + 



pd-l (1-TJd) ■ 



Remark 2 It is easy to see from Theorem U\ that 7Tq and vrg 1 ^^ are independent of the 
retrial processes including the two parameters 6 and d<i- Thus, the probing- server number 
will not have any influence on 7Tq and tt^^, but it can have influence on 7T^. and 
vrf ) /orfc> 1. 

In the remainder of this section, numerical examples are used to illustrate how the 
fixed point depends on the probing-server numbers d% and/or di- In the examples, we 
take A = 1, [i = 5 and 9 = 2. 
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When d 2 = 1, Table 1 shows how the fixed point depends on the probing-server number 



d\ > 2. It is easy to see from Table 1 that for a given d\ > 2, 7r[ and 7r^" ; are two 
monotone decreasing functions of > 1; and for a fixed k > 2, they decrease quickly to 
zero, as d\ increases. 



±w) 



Table 1: The fixed point is a function of d\ when d 2 = 1 





Aw) 

^0 


7T (/) 


(W) 
tt\ ' 


7T« 


Aw) 

n 2 


7T (/) 
^2 


di = 2 


0.1459 


0.8541 


0.0043 


0.0106 








dt = 5 


0.1110 


0.8890 














d x = 10 


0.0836 


0.9164 















When d\ = 1, Table 2 illustrates how 7Ti^ and vr^' 1 , for a fixed fc > 1, depend on 
the number d 2 > 2. It is seen that for a fixed k > 2, 7r[ and 7r£ are two monotone 
increasing functions of cfo > 2. Also, the crucial number Kd X) d 2 given in Proposition [2] is 
shown as follows: K\ 2 = ^15 = 0, i^i,8 = 1 and i^i,io = 5. Further, more numerical 
examples are used to illustrate that i^i,d 2 is a monotone increasing function of d 2 >2. 



Table 2: The fixed point is a function of d 2 when di = 1 







4" 


*™ 


4" 


AW) 
^2 


4" 


AW) 
^3 


4" 


d 2 = 2 


0.1667 


0.8333 


0.0911 


0.2887 


0.0484 


0.1509 


0.0276 


0.0898 


d 2 = 5 


0.1667 


0.8333 


0.1550 


0.6084 


0.1354 


0.5220 


0.1208 


0.4685 


d 2 = 8 


0.1667 


0.8333 


0.1799 


0.7330 


0.1717 


0.6786 


0.1626 


0.6413 


d 2 = 10 


0.1667 


0.8333 


0.1893 


0.7800 


0.1853 


0.7371 


0.1783 


0.7063 


d 2 = 15 


0.1667 


0.8333 


0.2028 


0.8473 


0.2045 


0.8197 


0.2006 


0.7983 


d 2 = 20 


0.1667 


0.8333 


0.2100 


0.8832 


0.2146 


0.8630 


0.2122 


0.8466 







An 


(W) 


4" 


AW) 


At) 






d 2 = 2 


0.0173 


0.0588 


0.0118 


0.0416 


0.0086 


0.0313 






d 2 = 5 


0.1106 


0.4323 


0.1034 


0.4062 


0.0980 


0.3864 






d 2 = 8 


0.1554 


0.6145 


0.1499 


0.5942 


0.1456 


0.5782 






d 2 = 10 


0.1724 


0.6836 


0.1677 


0.6663 


0.1640 


0.6524 






d 2 = 15 


0.1965 


0.7821 


0.1932 


0.7695 


0.1905 


0.7594 






d 2 = 20 


0.2093 


0.8340 


0.2067 


0.8242 


0.2046 


0.8162 
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Table 3 illustrates that for d± = 5 and a for a fixed k > 1, 7r^ and vr^-* are two 
monotone increasing functions of 62 > 2. Since the number d\ has a bigger impact on the 
fixed point than the number d2, it is seen that K^^ 2 = for d2 > 2. 



Table 3: The fixed point is a function of d 2 when d\ = 5 





^0 


4" 






ir (W) 






7T (/) 


d 2 = 5 


0.1110 


0.8890 





0.0967 














d 2 = 10 


0.1110 


0.8890 


0.0006 


0.3109 





0.0210 








d 2 = 15 


0.1110 


0.8890 


0.0041 


0.4589 





0.1456 





0.0209 


d 2 = 20 


0.1110 


0.8890 


0.0108 


0.5576 


0.0005 


0.3007 





0.1361 


d 2 = 25 


0.1110 


0.8890 


0.0193 


0.6267 


0.0029 


0.4297 


0.0004 


0.2899 


d 2 = 30 


0.1110 


0.8890 


0.0285 


0.6774 


0.0082 


0.5278 


0.0027 


0.4230 







7T (/) 


(W) 


7T (/) 




7T (/) 






d 2 = 5 
























d 2 = 10 
























d 2 = 15 





0.0008 


















d 2 = 20 





0.0498 





0.0140 





0.0029 






d 2 = 25 


0.0001 


0.1934 





0.1278 





0.0839 






d 2 = 30 


0.0010 


0.3483 


0.0004 


0.2941 


0.0002 


0.2539 







4 Expected Sojourn Time 



In this section, we use the fixed point to provide an iterative algorithm for computing 
the expected sojourn time E pd lj( 2 2 ] that a tagged arriving customer spends in the retrial 
supermarket model. 



TY_ L (k) 



and E 



If there are k customers in the orbit, then we denote by E 
the expected sojourn times that a tagged primary arriving customer spends in the re- 
trial supermarket model and that a tagged retrial arriving customer spends in the retrial 
supermarket model, respectively. At the same time, the probabilities that an arriving 
customer is primary and that an arriving customer is retrial are given by A/ (A + k9) and 
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kO / (A + k6), respectively. Thus we obtain 



k=0 



A 



where we take E 



X + k6 
0. 



:E 



r (P) 
~d 1 ,d 2 



(k) 



+ 



kO 



X + k9 



-E 



T 



(R) 



di,d2 



(k) 



(20) 



E 



(°). 

It is easy to see from Equation (|20p that we must first compute E 

(Ft) 1 

d2 (k) in order to compute E [?d 1)C £ 2 ]. To that end, we provide a useful relation 



- di ,d2 



(k) 



and 



between E 



T dud 2 (*) 



and £7 



r («) 



for k > 0. In what follows we consider two cases: 



E 



Case one: = 0. In this case, 

di 



r (P) 
L di,d 2 



(0) 



7T, 







7T 



(J) 



+ 



7T, 



(W0 



di 



7T 



E 



^il(l) • (21 



Note that 



7T, 



(0 



7T 



is the probability that a tagged primary arriving customer 



finds one idle server and no customer is in the orbit, thus the expected sojourn time of 



7T, 



(W0 



7T 



(W0 



is the 



the primary arriving customer is given by Similarly, 
probability that a tagged primary arriving customer finds that each server is busy and no 
customer is in the orbit. In this case, the tagged primary arriving customer has to enter 



the orbit so that this customer has the expected sojourn time E 



r (R) 

- di ,dl 



(1) 



Case two: k > 1. In this case, using a similar analysis to that in Case one we obtain 



E 



and 



E 



r (R) 

~di,d 2 



(k) 



7T 



(I) 



TT 



k+1 



+ 



TT 



(WO 



di 



7T 



(W) 
k+1 



E 



T Sd 2 (* + 1: 



(22) 



7T 



(I) 



d 2 



TT 



(0 
fc+i 



d 2 



+ 



7T 



(W) 



d 2 



TT 



(W) 
k+1 



E 



TZ, (k) • (23) 



Using Equations (|2ip to (|23p . the following theorem provides an iterative algorithm 
for computing the expected sojourn time. The proof is easy and is omitted here. 

Theorem 4 For k > 1 



E 



1 
/i 





[4 J) ] 


d 2 




Hi 


d 2 






(wo' 


d 2 




^k+l 
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E 



T T1 (o) 



i 

n 
+ 



7T, 



(I) 



TT 



(I) 



TT, 



(WO 



AW) 



d 2 



TT. 



(I) 



d 2 



/i 



1 - 



and 



E 



T T1 (*) 



1 



1 

+ - 



7T 



7T 



(W) 



TT 



(I) 
k+1 



di 



TT 



TT 



<h 



(W) 
k+1 



(W) 



TT. 



(W) 



d 2 



dl 



TT 



(I) 
k+1 



TT 



(I) 
k+2 



d 2 



TT 



(W) 
k+1 



d 2 



TT 



(W) 
k+2 



d-> 



Therefore, we have 

i — n v 



E 



T S 2 (*) 



fc=0 



kO „ 
H E 

x + ke 



L di,d 2 



TY\ (k) 



In what follows we consider a special case with d\ = 0I2 = d > 1. In this case, we have 



TT, 



(I) 



TT 



(I) 



(i - m) d - U 



and for k > 1 



for I > 



TT 



(I) 



TT 



(I) 
k+1 



7T 



7T 



l+l 



xd zd 
°l ~ d l+l- 



Thus we obtain that for A; > 1 



E 



T Z 2 (*) 



k9°k~l (k+l)0°k 



(k+l)6 l 



1 - 



k u k+l 



} ' 



rJS, (0) 



1 
/i 



and 



fc0 



+ 1)(9 



«52 +/> 



i-{5i-5 d 2 ) 

{ 5 k~ S k+l} { (fc+l)6^fc ~~ (fc+2)0^fc+l} 
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Therefore, we get 

E [ T d U d2] 



-E 

k=0 
1 

A* 



7# 



(1 - - 

oo 



d\,d2 



(k) 



+ 



k9 



X + k9 



E 



r (R) 

di,d2 



(k) 



+ P 



i-{sf-si} 



+ 



k=\ 



X P r 1 xd _ 1 ,d 

\ + ke x ko k - 1 (k + i)e k 



J k+l 



} {(k+l)9 5 k 



1- W.1-C2} 

fc6>/9 M^fe-l ~~ (fc+l)6»^fc 



; =i A + ^ 1-{8I-8I +1 } ■ 

In the rest of this section, we provide numerical examples to show how the expected 
sojourn time -E[Td li( 2 2 ] depends on the two probing-server numbers d\ and di- In these 
numerical examples, we take that A = 1, \i = 5 and 9 = 2. 

Figure 4 illustrates that the number d\ has a bigger influence on the expected sojourn 
time than the number c?2. It can be observe that d\ = 2 is a crucial point. If d\ > 
2, then influence of the number on the expected sojourn time is very weak because 
the three lines are almost identical. Figure 5 shows that the expected sojourn time is 
always monotone decreasing for c?2 > 1. It is seen that as the number d\ increases, the 
expected sojourn time has a diminutively decreasing rate for d2 < 100, but is can decrease 
quickly when > 1000. Figure 6 illustrates that the expected sojourn time is monotone 
decreasing for d\ = c?2 = d. It is seen from Figure 7 that the expected sojourn time is 
monotone decreasing for the two numbers d\ and cfe, and the number d\ plays a bigger 
role in the expected sojourn time than the number di- 

Remark 3 It is an interesting and difficult topic to study the role played by different 
probing-server numbers in performance analysis of supermarket models. Figures 4 to 7 
and Tables 1 to 3 for our numerical examples give insights into the role of two probing- 
server numbers used by the primary and the retrial customers, respectively. 



5 Exponential Convergence 

In this section, we provide an upper bound for the current location of the retrial super- 
market model. Based on this, we study exponential convergence of the current location 
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0.17 




Figure 4: The expected sojourn time depends on the number d\ 
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Figure 6: The expected sojourn time depends on the numbers d\ = d,2 
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d 

Figure 7: The expected sojourn time depends on the numbers d±,d; 
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to the fixed point. 

For the retrial supermarket model, the initial point S (0) can affect the current location 
S (t) for each t > 0, since the arrival, service and retrial processes in the retrial supermarket 
model are under a unified structure through a sample path comparison. To explain this, it 
is necessary to provide some notation for comparison of two vectors. Let a = (ai, 0,2,0,3, . . .) 
and b = (pi, 62, 03, . . .). We write a -< b if < for some k > 1 and ai < b\ for I ^ k,l > 1; 
and a ^ b if a& < for all fc > 1. 

Now, we can easily obtain the following proposition, while the proof is clear from a 
sample path analysis and is omitted here. 

Proposition 4 // S (0) * 5(0), then S (t) <S(t). 

Based on Proposition [H the following theorem shows that the fixed point ir is an upper 
bound of the current location 5 (t) for all t > 0. 

Theorem 5 For i/ie retrial supermarket model, if there exists some k such that Sjp (0) = 
or Sj[ (0) = ; i/ien i/ie current location S (t) for all t > /ias an upper bound vector 
ir, that is, S (t) ^ 7r /or a// 1 > 0. 

Proof: Let 

^ ) (o)=^,sf)(o)=.r ) , *>o. 

Then for each k > 0, 5^ (i) = 5^ (0) = tt^ and sf (i) = sf° (0) = ttJ^ for all 

■(I 
k 



t > 0, since S (0) = 7r is the fixed point for the retrial supermarket model. If si (0) = 



or (0) = for some k, then sjp (0) -« (0) or sjp (0) -< (0). Thus, we 



have S (0) ^ S (0). It is easy to see from Proposition |4] that si (i) ^ si (t) = 7r^ and 



S { ™ ] (t) 1 (t) = vrf for all fc > and t > 0, which lead to S (t) < vr for all t > 0. 
This completes the proof. ■ 
To show exponential convergence, we apply a potential function (or Lyapunov function) 
$ (t) which is defined as 



*(*) = £ - s« ( t)] + £ , fc - ^ (*)] , (24) 

k=0 k=0 

where {w^} and {ufc} are two underdetermined constant sequences with wo = 1. 

The following theorem measures the distance of the current location S (t) for t > to 
the fixed point 7r, and illustrates that the L\ distance to the fixed point from the current 
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location is very close to zero with exponential convergence. Hence, it shows that from any 
suitable starting point, the retrial supermarket model can quickly converge to the fixed 
point, that is, there always exists a fixed point in the retrial supermarket model. 

Theorem 6 For t > 0, 

$(t)<c e~ 5t , 

where cq and 5 are two positive constants. In this case, the potential function (t) is 
exponentially convergent. 



Proof: It is seen from (|2~1| that 

oo 

*(*)=!>* K" - 4" «]+!> K> - sP w] . 

which leads to 



k=0 



k=0 



fc=0 



fc=0 



It follows from Equations (|T|) to (jlOp that 



fc=o 



d 2 



- E ^^ 5 i 7) (*) = E Xw * \ s k ] (*)1 1 + E ^ N 7) 2 - E (*) 



fc=0 



k=l 



k=0 



and 



fc=o 

Thus, we obtain 
d 



E -4^ } (t) = - E ^ Ufc (*)] 1 - E H T) (*) 



d 2 



fc=0 

oo 



fc=l 



+ e a k - , fe+1 ) (t)] di + x: /^r w , 



fc=0 



k=0 



— *(t) = ^ A (ufc - u fc 

fc=0 

oo 

+ E A ( Ufc ~ Vfe+i 



fc=0 



fc=i 

, oo 



k=0 



Let 



H 7) 




di 


= c k (t) 








d 2 


= 4 (t) 
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d 

— $ (t) = ^ A K - Ufc ) c fe (t) [vr fc - 4 7) (t) 

fc=0 

oo 

+ M (w k - 4 (t) [vr fe - 4 7) (*) 
fc=i 

oo 

+ ^2^(vk~ V k + l) 9k (t) TTk ~ (t) 

k=0 

oo 

+ ^2^( v k~ w k ) h k (t) vr fc - S { k W) (t) 

k=0 

=A (W - V ) C (t) 7T - Sjf* (t) 
oo 

fc=l 

oo 

+ t A ~ v k+i)9k (t) + fJ,(v k - w k ) h k (t)] ir k - S k W) (t) 



and 



Then 



Let 



for k > 1 



and for / > 



fc=o 

u> = 1, 
A (w Q - v ) c (t) < -Svo, 

A (w k - v k ) c k (t) + kO (w k - v k _i) d k (t) < -Sv k ; 

A (yi - vi +1 ) gi (t) +n(vi- w t ) h t (t) < -Swi. 



Then 



7 I oo oo 

±» w < -t E «* K» - 4" w] + E <* K"' - w 



.fc=0 



fc=0 



This gives 



thus we obtain 



d_ 

lit 



$ (t) < (t) , 
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This completes the proof. 



Remark 4 We provide an algorithm for computing the two underdetermined constant 
sequences {uik} and {v^} with Wq = 1 as follows: 
Step one: 

w = 1. 



Step two: 
Step three: 



A 



and 



Ac (t)-<T 

5 + Xv go (t) + fj, (vp - 1) h (t) 
Xgo (t) 

Auici (i) + 0UQ(ii (t) — (5^1 



Aci (t) + Odi (t) 
Step four: for k > 2 

= Swk-x + Xvk-igk-i (t) + At Q fc _i - w fc _i) fefc-i ft) 

and 

_ A-u fc Cfc (t) + kOvk^idk (t) - 5v k 
Wk ~ Xc k (t) + Atfdfc (t) 

6 A Lipschitz Condition 

In this section, we apply the Kurtz Theorem to study the density-dependent jump Markov 
process corresponding to the retrial supermarket model with Poisson arrivals, exponential 
service times and exponential retrial times and with two different probing-server numbers. 
This will lead to the Lipschitz condition under which the fraction measure of the retrial 
supermarket model weakly converges to the system of differential equations. 

The retrial supermarket model can be analyzed by a density-dependent jump Markov 
process, which is a Markov process with the population size n. Kurtz's work provides 
a basis for density-dependent jump Markov processes in order to relate infinite-size sys- 
tems of differential equations to corresponding finite-size systems of differential equations. 
Readers may refer to Kurtz [16] for more details. 

In the retrial supermarket model with n servers, we denote by nl the numbers of 
busy servers with at least k customers in the orbit, and nl the numbers of idle servers 
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with at least k customers in the orbit. Hence < n^\n^ < n for k > 0. Note that 



and n^' are independent of time t > 0, thus they are different from n£ v 1 (t) and 
n^P (t) defined in Section 2 for the retrial supermarket models. 
We write 

= {(w, k) : < k < n} 

and 

EW ={(I,k):0<k<n}, 

In the state space E n , the density-dependent jump Markov process {X n (t) : t > 0}, defined 
in Section 2, for the retrial supermarket model contains four classes of state transitions: 



(W) 



• ;ii > „ nm. ^ W ,i)^(W,i + l), 
Chx^ two: E y n ' s^E K n>: (W,i) -»• (J,i), 



Class one: £^ ; ^ ^ 



Class three: E, 



V) 



Class four: r. E^: 



(W,i), and 
(W,i-1). 

Note that the state transitions a^, and express a primary arrival, a service completion 
and a retrial arrival, respectively. 
For k > we write 



and 



Let 



4 7) (n) = S 



i W) (n) = £ 



r?. 



U) 



n 



5*0 (n) 



n 



(n) 



S fe (n) = ( S (n) , a? (n)), 



and for A; > 1 



S(n) = (S (n),S 1 (n) , S 2 (n) , . . .) . 

Note that the states of the density-dependent jump Markov process {X n (t) : t > 0} can be 
normalized and be interpreted as measuring population densities: S (n) = (So (n) , Si (n) , S 2 (n) , . . .). 
The transition rates of this Markov process also depend on these population densities. 
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(n) 

1(W,k)^(I,k) ~ n l- 
(n) 



(n) 

l(I,k)->-(W,> 



Let {x„ (t) : t > 0} be a density-dependent jump Markov process on the state space 
.En whose transition rates corresponding to the above four cases of state transitions are 
given by 

1 ( (w,k)^(W,k+l) = n h _a fc+i (4 W) ( n )) = nXs k W) ( n ) » ( 25 ) 
h^k{s^ ) {n))=n^\n), (26) 

n/3 fc a fc (4 7) (n)) = nAs£° (n) , (27) 

= n/3 fc ^ fc _! (4 7) (")) = ( n ) • ( 28 ) 

In the retrial supermarket model, X n (t) is an unsealed process which records either the 
number of busy servers or the number of idle servers with at least k customers in the orbit 
for k > 0. Using Chapter 7 in Kurtz [16] or Subsection 3.4.1 in Mitzenmacher [28J, the 
Markov process |x„ (t) : t > o| with the transition rates, described in Equations (|25|) to 
f[28|) . is given by 

X n (t) = X n (0) + £ / b n jf ' A d ^ , (29) 

where Y b (x) are independent standard Poisson processes and l b is a positive integer with 
h < 5? < +00 for € E, where 

E = {(W,ife) _^ (W, fc + 1) , (W, fc) 4 (I,k),(I,k) (W, k) ;(I,k + 1) (W, k) : k > 0}. 

Clearly, the state of the density-dependent jump Markov process in Equation f|29|) at time 
i is determined by the starting point and the transition rates in Equations (|25p to (|28p . 
which are integrated over its history. 
Let 

F(y)= Y, hPbiv), (30) 

beE(y) 

where 

E (y) = {b € E : A state transition 6 is directly related to state y} . 

That is, the state transition b from state y contain two classes: 1) y and y for 
entering state y, and 2) by and y b^ for leaving from state y. Taking X n (t) = n~~ 1 X n (t) 
which is an appropriate scaled process, we have 

"t \ ft 



X n (t) = X n (0) + V l b n~ l Y b (n [ b (X n («)) du) + [ F (X n («)) du, 
b€E V Jo J Jo 



(31) 
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where Yj> (y) = Yj, (y) — y is a Poisson process centered at its expectation. 
Let X (t) = linin._j.oo X n (t) and xq = lim^^eo X n (0), we obtain 



X(t)=x + f F(X(u))du, t>0, 
Jo 



(32) 



due to the fact that 



lim -% ( n I ft (X n (_)) du 

n— too fi 







by means of the law of large numbers. In the retrial supermarket model, the deterministic 
and continuous process {X (t) ,t > 0} is described by the infinite-size system of differential 
equations (|7|) to (fTUj) . or simply, 



jX{t) = F{X{t)) 



(33) 



with the initial condition 

X (0) = xq. (34) 

In the rest of this section, it is necessary to express the function F (y) in order to set 
up a Lipschitz condition. To that end, using the law of large numbers we have 

Tik = lim Sk (n) , k > 0. 

n— too 

Let 



n 



{^vsvr^vr;...}. 



'0 ' "1 ' "0 ' 

Then it is easy to see from the system of differential equations ([7]) to (IIPD that the function 
F (y) is given by 



F{y) 



7T, 



(I) 



A 



7T, 



7T 



(/)' 



di 



(W) 







7T 



A 
+A 



7T 



(W) 

fc-1 
(0 



7T, 



(W) 



4 



7T 



d 2 

+ 

di 



(WO 



•r (0 

if _ = TTo > 

if y = 7r£ for k > 1, 



7T 



d 2 



(wo 



+ + 



7T 



(0 
fc+1 



<«2 



7T 



(wo 



if y = vr^ 



if y = 7r£ for A; > 1. 



(35) 



Now, we consider uniqueness of the limiting deterministic process {X (t) ,t > 0} with 
and (|34|) . or uniqueness of solution to the infinite-size system of differential equations 
([7]) to (|10p . To that end, a sufficient condition for uniqueness of solution is Lipschitz, that 
is, there exists a number M > such that for y, z £ fi 



|F(y) - F(z) I < M||y - z\\. 
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In general, the Lipschitz condition is standard and sufficient for uniqueness of solution to 
the finite-size system of differential equations; while for the countable infinite-size case, 
readers may refer to Theorem 3.2 in Deimling [7] and Subsection 3.4.1 in Mitzenmacher 
[28] for some useful generalization of the Lipschitz condition. 

The following theorem shows that the retrial supermarket model satisfies the Lipschitz 
condition for analyzing uniqueness of solution to the infinite-size system of differential 
equations (J7]) to (fTUI) . 

Theorem 7 The Lipschitz condition is satisfied in the retrial supermarket model with 
Poisson arrivals, exponential service times and exponential retrial times and with two 
different probing- server numbers. 



Proof For the retrial supermarket model, it is easy to see from (I35p that we need 
to consider eight types of pairs (y, z) in order to check that there exists a bigger M > 
such that \F (y) — F (z) \ < M\\y — z\ \ for any y, z E fi. In what follows we only prove two 
types of pairs (y, z), while all the others can be proved similarly. 



Let y 



4 and z 



7T^, ' for k > 1. Then 



\F(y)-F(z)\<\\ 



7T, 



7T 



(0 



I | I (W) 
l+MFo 



k6 



7T 



dz 



Note that 



which leads to 



7T 



(0 



—S dl 
k6 k ~ l 



d 2 



kO 



7T 



(I) 



d 2 



Since 4^ > 4^ and 4^ > 4 , we have 



\F(y)-F(z)\<\ 

< A 

< A 



7T, 



7T, 



(0 



(0 



7T 



(0 



\+v {4 



(W) 



7T, 



(0 



+ K + A 

+ v.*™ + A 



7T 



7T 



7T, 



(WO 
fe-i 

(wo 



due to the fact that ir^ > tt { ™\ > vrj^ for k > 2. Since vr^ > 4-1 > 4° for k ^ 2 ' 



we have 



7T, 



(0 



4" < 4 J ) - 



k ' 
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which leads to 



A 



7T, 



(I) 



dl 



7T, 



(WO 



di 



« J 1 ) 



> 



(I) 



di 



7T, 



(W) 



dl 



Taking 



we obtain 



■7T n - 7T 



Mi = 



7T (/) - 7T (/) 



+ /U7TQ + A 



7T, 



(wo 



di 



\F(y)-F(z)\<M 1 



Mi|y-z|. 



Let y = 71"! and z = iv k for A; > 2. Then 



|F(y)-F(z)|<A| 



r (W0 



+ |26> 
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di 
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Note that for At > 2 
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which follows 
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3A 



Mo 
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7T 
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|F(y)-F(z)| <M 2 



7T 



(W) _ (W) 



M 2 \y-z\. 



Similar to the above analysis, for the ith type of pairs (y,z) there exists a number 
Mi > such that 

|F(y)-F(z)|<M||y-z||. 
Taking M = max {Mj : i = 1, 2, . . . , 8}, we obtain that for any y, z G fl, 

\F(y)-F(z)\<M\\y-z\\. 
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This completes the proof. ■ 
Based on Theorem the following theorem directly follows from Theorem 3.13 in 
Mitzenmacher 128 . 



Theorem 8 In the retrial supermarket model with Poisson arrivals, exponential service 
times and exponential retrial times and with two different probing- server numbers, we have 

lim sup \X n (u) — X (u) I = 0, a.s., 

n ^°° u<t 

where {X n (t)} and {X (t)} are given by \31\) and A32\), respectively. 

Proof For the retrial supermarket model, Theorem [7] shows that the function F (y) 
for i/SS] satisfies the Lipschitz condition. Based on this, it is easy to take a subset O* C Q 
such that 

{X (u):u<t}c n* 

and 

sup F (y) < +oo. 
yen* 

Thus, this proof can be further completed by means of Theorem 3.13 in Mitzenmacher 
|28j. This completes the proof. ■ 
Using Theorem 3.11 in Mitzenmacher [28] and Theorem [HJ the following theorem pro- 
vides an upper bound over the time interval [0, t] for the expected sojourn time that an 
arriving tagged customer spends in an initially empty retrial supermarket model. 

Theorem 9 For the retrial supermarket model with Poisson arrivals, exponential service 
times and exponential retrial times and with two different probing- server numbers, the 
expected sojourn time that a tagged arriving customer spends in an initially empty system 
over the time interval [0, t] is bounded above by 

k=o 

where o(l) is understood from the population size n — > oo. 

7 Concluding remarks 

In this paper, we have presented a novel approach for studying super-exponential solu- 
tion of the retrial supermarket model with Poisson arrivals, exponential service times and 



rf ) (k) 



+ o(l), 
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exponential retrial times, and with two different probing-server numbers. We describe 
the retrial supermarket model as a system of differential equations by means of density- 
dependent jump Markov processes, and obtain an iterative algorithm for computing the 
fixed point of the system of differential equations. Based on the fixed point, we analyze 
the expected sojourn time that a tagged arriving customer spends in this system, and use 
numerical examples to indicate different influence of the two probing-server numbers on 
system performance including the expected sojourn time. Furthermore, we analyze expo- 
nential convergence of the current location of the retrial supermarket model to the fixed 
point, and apply the Kurtz Theorem to study density-dependent jump Markov process 
given in the retrial supermarket model. This leads to a Lipschitz condition under which 
the fraction measure of the retrial supermarket model weakly converges to the system of 
differential equations. This paper gives a new understanding of how the workload probing 
can help in load balancing jobs in retrial supermarket models. 

Our approach given in this paper is useful in the study of retrial supermarket models 
with applications to, such as data centers, multi-core servers systems and cloud compu- 
tational modeling. We expect that this approach will be applicable to the study of more 
general retrial supermarket models, such as non-Poisson arrival processes, non-exponential 
service times, and retrial supermarket networks (for example, each server contains c iden- 
tical processors). It will also open a new avenue to quantitative evaluation of more general 
supermarket models and analysis of retrial queueing networks. 
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